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Structure of confining string in terms of the topological charge density and the action density is 
studied in 5*17(2) Yang-Mills theory on the lattice using HP u-model embedding approach. We 
find that the confining flux tube noticeably suppresses both the topological charge and the action 
densities. Beyond the string formation length the string cross section in terms of these quantities 
is well described by a Gaussian profile. In both cases the squared string width is found to be a 
logarithmic function of the string length confirming the Luscher widening of the chromoelectric 
string. Characteristic string scales in terms of the topological and action densities are estimated as 
well. 
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I. INTRODUCTION 

It is commonly believed that the phenomenon of the 
quark confinement in QCD happens due to formation 
of the chromoelectric string spanned between quark and 
anti-quark. In fact, various numerical simulations in 
QCD both in quenched and in unquenched limits indicate 
that in the QCD vacuum the fundamental color charges 
(quarks) are connected by a string-like structure [l| (see, 
e.g., Ref and references therein), which can be seen 
directly in gluonic energy and action density [H, Q| as well 
as in other related quantities [5] . Physically, the string is 
formed by squeezed fluxes of chromoelectric fields com- 
ing out of quarks. The chromoelectric string possesses a 
non-zero tension providing steady attracting force which 
makes the quarks confined into colorless hadrons. 

There are various effects caused by the strong chro- 
moelectric fields of the QCD string. Besides the en- 
hancement of the gluonic energy density, the string is 
characterized by the suppression of the gluonic action 
density @, as well as by the suppression Q of the (A 2 )- 
condensate Q. Some features of the QCD string are con- 
sistent with the dual superconductor model of the color 
confinement @ which attributes the color flux squeezing 
to a (dual) Meissner effect originating from condensa- 
tion of particular ( "monopole-like" ) gluon field configura- 
tions [10( 1 . Naturally, chromoelectric string formation has 
a great impact on both gluons and quarks, however, in 
latter case the confinement of color manifests itself in yet 
another way. There are numerical indications that the 
chiral symmetry, originally broken by non-perturbativc 
effects in the low temperature QCD, is (partially) re- 
stored inside the string [ll[ ■ The effect is revealed by lo- 
cal suppression of the chiral condensate (ip^) in the vicin- 
ity of the string. The same effect is also observed numeri- 
cally in the neighborhood of individual heavy quarks [r3 | , 
indicating that the strong chromoelectric fields tend to 
restore the chiral symmetry. Since the chiral proper- 
ties of the vacuum are tightly related to its topology, 
one expects generically that strong chromoelectric field 
of the string severely affects the topological charge den- 
sity. Indeed, the anti-correlation between the QCD string 



and the topological charge density was discovered in nu- 
merical simulation of Ref. [l3| in the so called "cooled" 
SU(3) vacuum. The investigation was developed further 
in Ref. [14j with the help of the same cooling method 
which is a gradual transformation of the quantum ( "hot" ) 
fields towards the classical ( "cold" ) vacuum. The cooling 
is aimed to reduce the topologically trivial effects of di- 
vergent zero-point fluctuations. This important result - 
obtained with the simplest lattice definitions of the topo- 
logical charge - is a first numerical indication that the 
QCD string suppresses topological fluctuations. In this 
paper we investigate, in particular, the structure of con- 
fining string in terms of the topological charge density. 
The problem is conceptually and technically challenging 
and therefore it calls for a few physically motivated sim- 
plifications. First, we refrain from using the cooling ap- 
proach because it certainly wipes out local topology of 
the real (uncooled) vacuum. Indeed, the fluctuations of 
the topological charge inside the chromoelectric string 
were found to depend crucially upon the cooling param- 
eters [l3|, [HI and therefore are affected by uncontrollable 
systematic errors. Other essential simplifications used in 
this paper are as follows: 

Quenched vacuum. For our purposes it is natural 
to consider quenched approximation because dynamical 
quarks screen the test color charges and break the string. 
As a result the correlation between the string and any lo- 
cally defined quantity disappears at sufficiently large sep- 
arations between the test charges. Thus the choice of the 
quenched vacuum is crucial, as it allows to discriminate 
between the string breaking and the actual correlation 
effects. 

Reduced number of colors. It is advantageous to uti- 
lize the known similarity between the SU (2) and SU (3) 
gauge models. Since nonperturbative aspects are alike in 
these theories, we naturally expect that the behavior of 
the topological charge density around the string in the 
realistic SU(3) case is qualitatively the same as for the 
SU(2) gauge group studied in our paper. Apart from 
the obvious computational simplifications, restriction to 
SU(2) gauge theory allows us to address unambiguously 
the next urgent problem of 
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Ultraviolet (UV) fluctuations. Indeed, simply on di- 
mensional grounds one expects that the vacuum topo- 
logical density is power-like divergent (q 2 ) ~ 1/a 8 , a be- 
ing the lattice spacing (the inverse UV cutoff). How- 
ever, the leading UV divergence should not be sensitive 
to the quark-antiquark separation. Note that this is not 
so in the case of subleading power corrections, which de- 
pend also upon the infrared (IR) scale (denoted Aqcd 
generically) . Therefore in order to study the correla- 
tion of the string and the topological density it would be 
advantageous to remove the leading (and only leading) 
power divergence in (q 2 ) since it gives identically vanish- 
ing but otherwise numerically destructive contribution to 
the correlation function. 

In view of the above points the choice of the topological 
density operator, which is free from additional parame- 
ters, becomes, in fact, almost unique. Indeed, the cool- 
ing approach was rejected for reasons explained above. 
As far as the fermionic definition [l5| is concerned, in 
the numerical implementations it is either not parameter 
free or is not capable to filter out exclusively the lead- 
ing power divergence. Indeed, the essence of the con- 
struction is to sum up the individual contributions of 
eigenmodes of overlap Dirac operator, ranked with cor- 
responding eigenvalues. It is known [l6| that the full 
tower of eigenmodes indeed reproduces ~ I/a 8 behav- 
ior, while the restriction to lowest eigenvalues [TtJ intro- 
duces explicit and a priori ambiguous cut on the leading 
singularity. Therefore the essentially unique alternative 
construction which fits the above requirements is HP 1 
cr-model embedding approach [H, which is to be re- 
viewed in Section [TTJ Here we only mention that this 
method indeed provides parameter free definition of the 
topological density in which the leading power divergence 
is removed. Note that the removal of leading perturbative 
contributions from various observables in HP 1 embed- 
ding approach had not been proved rigorously, however, 
in all examples considered so far it turns out to be the 
case. Moreover, the subleading power corrections seem 
to be left intact by the corresponding gauge covariant 
HP^projection procedure. What is even more impor- 
tant is that near the continuum limit the HP 1 -projected 
fields capture the majority of non-perturbative aspects of 
SU(2) Yang-Mills theory and reproduce, for instance, the 
full tension of the confining string and (quenched) quark 
condensate (see Refs. flU TTolj for further details). Be- 
ing numerically superior, the HP 1 embedding approach 
is naturally our method of choice to investigate the topo- 
logical aspects of the confining string and to study its 
geometrical properties. For discussion of other methods 
which allow to "filter out" topological properties of the 
vacuum out of the ultraviolet noise see Ref. [20| . 

However, the advantages of our method are not coming 
for free, the price is the intrinsic non-locality: the HP 1 
cr-model fields constructed from given SU(2) gauge back- 
ground are not local in terms of the original gauge poten- 
tials. However the characteristic range of non-locality is 
certainly less than < 0.3 fm and is comparable with the 



expected [l| half- width of the flux tube as it is determined 
by the string action density profile. We conclude there- 
fore that the non-locality of our method should not be 
essential for large enough quark-antiquark separations. 

The knowledge of the string profiles immediately re- 
veals an important feature of the chromoelectric string 
known as the string widening. This rather general phe- 
nomenon occurs in various gauge and spin systems pos- 
sessing the string-like excitations. As it was predicted 
by Luscher plj , quantum fluctuations force the string to 
fluctuate transversely in such a way that its geometrical 
center follows the Gaussian distribution. The squared 
width of the Gaussian distribution (string width) is ex- 
pected to be the logarithmic function of the string length, 
S 2 (R) oc In R, which in quenched QCD is the distance be- 
tween heavy quark-antiquark pair. 

The logarithmic widening of the string was accurately 
seen [22[ in the %2 gauge model in (2 + 1) dimensions, 
and the signatures of the widening are reflected in the 
form of the potential between test sources connected by 
the string [23|]. For other gauge theories the results are 
still inconclusive. For instance, the data for the string 
in gluodynamics width obtained via the energy density 
profiles [13 is consistent both with the constant behavior, 
S(R) » const, and with the logarithmic scaling for R > 
1 fm. 

In the numerical simulations the string profile is usu- 
ally overshadowed by the ultraviolet noise, which may be 
reduced within the certain models of the vacuum struc- 
ture. For instance, the dual superconductor scenario of 
quark confinement Q suggests that the confining string 
formation is governed by the dynamics of Abelian fields 
calculated in a certain Abelian gauge. While this ap- 
proach consistently describes particular features of the 
chromoelectric string, the string width is found to be 
length-independent both in quenched Q and full [24[ 
QCD, thus contradicting the logarithmic scaling law. 
Another popular approach [25j assumes that basic non- 
perturbative properties of gluons are encoded in the cen- 
ter degrees of freedom , th e corresponding string profile 
was calculated in Ref. [2g|. It was found that the string 
width itself (not its square) is a logarithmic function 
of the quark-antiquark distance, again contradicting the 
Luscher law. Therefore the results available in the lit- 
erature indicate that neither Abelian nor center degrees 
of freedom preserve the transverse non-perturbative fluc- 
tuations of the confining string, while the simulations 
within the full theory give inconclusive results. In this 
paper we investigate the string broadening using HP 1 - 
projected fields, which allow us to obtain a clean picture 
of the confining string both in terms of the topological 
density and in terms of the gluonic action distribution. 

The structure of the paper is as follows. The essentials 
of HP 1 cr-model embedding approach are reviewed in Sec- 
tions |TT1 while in Section IIIII we discuss the relevant ob- 
servables and theoretical expectations. In Section ITVl we 
describe the technical details of our measurements. Then 
the profiles of the QCD string in terms of the topological 
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and HP 1 -projected action densities are presented in Sec- 
tions [VjVll respectively, where we also discuss the degree 
of the string broadening at various quark- antiquark sep- 
arations. Section IVlII is devoted to the discussion of our 
results. 



II. REVIEW OF HP 1 EMBEDDING 

In this Section we briefly review the HP^ er-model em- 
bedding approach using continuum notation, for precise 
lattice definitions and further details see Refs. [la. fTsj] . 
The most straightforward way to introduce the idea of 
the method is to consider first Yang-Mills theory with 
gauge group SU(2) (group of unit quaternions) in the 
limit, in which the gauge potentials are smooth func- 
tions of space-time coordinates. Then it is known (2?| 
(see, e.g., Ref. [28| for extended discussions) that A^ can 
be represented as 

N 

= (q\dfj.\q) = ^Qid^ , (1) 
i=0 

where \q) is (N + 1) x 1 normalized (q\q) — 1 col- 
umn vector with quaternion valued entries g,, bar de- 
notes quaternionic conjugation and N is finite integer 
number, the upper bound on which is unimportant for 
the present discussion. Gauge covariance of (fT]) implies 
that \ q) = [qa,...,q N ] T and \q)v= [q v, c/atu] t , where 
v is unit quaternion, vv = 1, are physically equivalent 
and therefore the set of gauge inequivalent fields | q ) de- 
scribe the quaternionic projective space HP^, convenient 
parametrization of which is provided by gauge invariant 
(JV+1) x (JV+1) projection matrices P = \q)(q\. There- 
fore in the limit of smooth fields the SU(2) Yang-Mills 
theory could be reformulated in terms of the non-linear 
cr-model with target space RP N . Note that in the partic- 
ular case N = 1 target space geometry simplifies greatly, 
HP 1 = S* 4 , and one has 

P=\q){q\ = \{l + 1 A n A ), A=l,...,5, (2) 

where n A is unit five-dimensional vector and j A are the 
conventional Dirac matrices. We use the same symbol 
q to denote both the quaternions and the topological 
charge density. We do not expect a confusion with nota- 
tion since the quaternions are used only in this Section 
and they are always associated with the bra-ket symbols. 

The essence of HP^ cr-model embedding approach is 
to consider Eq. (fT]) within the Euclidean formulation of 
quantum SU (2) Yang-Mills theory and to quantify its ap- 
proximate validity for given gauge background. Namely, 
for any fixed gauge configuration A^ one finds the "near- 
est" HP^ a- model fields {|c/ x )}) which minimize the 
functional norm 

\\A-{q\d\ q )\\=F(A,\ q ))<F(A,\q')), (3) 



In fact, the problem ([3]) is well suited for the lattice stud- 
ies and it turns out that the minimizing configuration 
is essentially unique, at least for the gauge po- 
tentials numerically generated on the lattice near the 
continuum limit. Note that the only free parameter 
at this stage is the rank N of the cr-model considered. 
Naively, one could expect that the quality of approxima- 
tion A sa ( q \d\ q) quantified by F(A, | q )) varies strongly 
with N, however, this turns out not to be the case. To 
the contrary, the functional distance F(A, | q)) appears 
to be independent on the er-model rank and hence there 
are no practical reasons to go beyond the simplest N = 1 
case. In turn the restriction to N = 1 allows to give 
natural and computationally superior definition of the 
topological charge density, which undoubtedly integrates 
to an integer number. 

At this point it is natural to introduce the notion of 
HP 1 -projection. Namely, one could replace the origi- 
nal potentials with corresponding cr-model induced gauge 
connection 

A, Af> = (q\d„\q), (4) 

and then compare the physical content of the original and 
projected configurations. It turns out that A^ F fields 
are extremely weak, the corresponding curvature is by 
order of magnitude smaller than that of A^ configura- 
tions. However, the detailed analysis show that (we are 
necessarily brief here, see Refs. [3, [l9| for details): 

• projected theory is still confining, moreover, the 
projected string tension cthp exactly reproduces the 
full one in the continuum limit; 

• topological aspects of the original fields are pre- 
served. This could be verified by various methods 
including modern fermionic ones. 

• investigation of the local observables, like topologi- 
cal charge or action densities, suggests that only the 
leading perturbative power divergences are washed 
out by HP^projection. In particular, the method 
allows to estimate rather precisely both the gluon 
condensate and its leading (quadratic) power cor- 
rection, the magnitudes of which are in full agree- 
ment with the existent literature. Moreover, it was 
found that the topological fluctuations behave quite 
non-trivially in the continuum forming the perco- 
lating three-dimensional domains of sign coherent 
topological charge in the limit of vanishing lattice 
spacing. The latter observation is again in full 
agreement with modern literature [H, IT?], H(| • 

The overall likely conclusion is that the HP 1 embed- 
ding method is the perfectly gauge invariant way to filter 
out only the leading UV divergent part from various ob- 
servables leaving intact the non-perturbative content of 
the original gauge configurations. However, there is a 
particular disadvantage: the cr-model fields constructed 
this way are non-local in term of the original gauge poten- 
tials. At present, there are two possibilities to estimate 



4 



the range of non-locality. First, the heavy quark poten- 
tial extracted from HP 1 -projected Wilson loops (Figure 
9 in the second paper of Ref. [Hi]) appears to be concave 
at distances \x\ < 4 (lattice units) for gauge configura- 
tions thermalized at bare coupling (3 = 2.60. Secondly, 
the HP 1 -based topological density correlation function 
(lolx) reveals a positive core at small distances, the size 
of which scales in physical units [l9( . Both these observa- 
tions indicate that the characteristic range of non-locality 
is (see Section HVf 



R 



non— local 



V^hp < 0.5 . 



(5) 



Note that although this number is parametrically large, 
Eq.© provides only an upper bound. We expect that 
for large enough Wilson loops the corrections due to the 
non-locality are negligible provided that the linear extent 
and the width of the string is greater than R no n— local' 
For this reason in the numerical investigations below we 
systematically exclude data points for which the relevant 
geometrical characteristics are smaller than i? n on-iocai- 



III. OBSERVABLES AND THEORETICAL 
EXPECTATIONS 

Before going into details of the numerical simulations 
let us discuss the observables to be measured. As far 
as the topological density (q 2 ) and HP 1 -projected action 
density (s) are concerned, their construction is fully de- 
scribed in Refs. [H, Eil and will not be repeated here. 
The prime objects of our study are the dimensionless 
quantities 



Qe(x) 



Sc(x) 



[W(C)q 2 (x) 
(W(C)) (q 2 ) 



[W(C)s(x) 
(W(C)) (s) 



1, 



(6) 



(7) 



which are the expectation values of squared topological 
density (q 2 ) and HP 1 -projected action density (s) in the 
presence of the external color sources represented by the 
rectangular C = T x R Wilson loop W(C) in the funda- 
mental representation. Note that W(C) is to be evaluated 
also on HP 1 -projected potentials 



W(C) = Tr Pexpjy dtx^A™} 



(8) 



where t is natural parametrization of C and dot denotes 
derivative with respect to t. The normalization of both 
quantities Qc{x) and Sc(x) is such that in the case of neg- 
ligible correlation with the confining string both quanti- 
ties vanish. Furthermore, suppression or enhancement 
of the topological/action density fluctuations makes the 
corresponding quantity negative or positive, respectively. 




FIG. 1: (Color online) The geometry of our simulations. The 
topological/action densities are calculated at the point (r,h), 
where r and h are, respectively, the transverse and longitu- 
dinal coordinates with respect to the geometrical center C of 
the rectangular T x R Wilson loop. 



It is amusing to note that the HP 1 -projected Wilson 
loop has a pure geometrical interpretation. Indeed, using 
Eqs. (|2I4|) one can show that 

W(C) = Tr Pexpj § J c dt cj ab lo ab (t) } = (9) 
det^d AB d t - ^wab} 



1/2 



where a AB = f[7 A ,7 B ], A,B= 1,...,5 and 



,AB 



\{n A h B 



n B h A ) 



is the instantaneous rota- 



tion speed of the color vector n A {t) naturally assigned to 
the contour C (see, e.g., Ref. HU f° r detailed derivation). 
Therefore the HP^projected Wilson loop is nothing but 
the spin factor associated with spinning particle propa- 
gating in fictitious five-dimensional color space, n A being 
the normalized tangent to its trajectory. The spinning 
particle interpretation of the Wilson loop might have a 
deep consequences, however, what precludes further anal- 
ysis is the unknown measure of integration over n A field. 
Indeed, at present we only know how to generate n A con- 
figurations numerically and have no idea what is the key 
feature of the measure T>n, which provides an area law 
for (W(C)}. In principle, one could speculate that the 
measure T>n might be taken standard and the area law 
would still be captured provided that the flat metric S AB 
in © is replaced by some effective one g AB describing 
geometry of ambient 5-dimcnsional space. This way of 
thinking resembles somewhat the AdS/QCD framework 
(sec, e.g., Ref. 32]), however, it would lead us too away 
from the purposes of present discussion. 

The geometry relevant to the correlation functions ((6|) , 
([7]) is shown in Figure [TJ For the rectangular T x R Wil- 
son loop the chromoelectric string appears between the 
static quark- antiquark pair at the moment t = 0. Once 
created, the string develops in full at the time t = T/2 
when the numerical evaluation of the string profile is per- 
formed. Finally, the string disappears at the time t = T 
when the quark and the antiquark are annihilated. For 
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large enough T the problem evidently possesses spacial 
cylindrical symmetry, which suggests to use coordinates 
x = (r, h), where r and h are, respectively, the transverse 
and longitudinal distances to the geometrical string cen- 
ter C. 

As is well known, at the quantum-mechanical level 
the string must fluctuate around classical (minimum ac- 
tion) configuration. Liischer [2l[ has obtained long ago 
a quantitative estimation of the quantum fluctuations of 
a bosonic string. According to Ref. [21] one generally 
expects that the deviation of the string position with 
respect to its classical limit should follow a Gaussian dis- 
tribution with a calculable width. The quantum fluctu- 
ations of the string should also affect distributions of all 
other quantities associated with the string. In particular, 
the profiles of the action and topological densities should 
become noticeably wider due to the string fluctuations 
compared to the that for the classical (non-fluctuating) 
string. Therefore, one expects generically that the trans- 
verse slices of the correlation functions ([6]), (O at the 
geometrical center ft, = of the string should follow the 
Liischer scaling law: 



X c (r,R) 



-° x{R) 6XP {-^)}' 



(10) 



where X stands for either the topological charge or the 
action, X = Q, S . This distribution is characterized by 
two quantities which are the amplitude Cx and the width 
Sx of the transverse string fluctuations. Note that both 
these quantities depend in general on the length of the 
string R. In particular, the squared width of the bosonic 
string is known [2l| to increase logarithmically for suffi- 
ciently long strings: 



d 2 



R 



DC ore 

H x 



(11) 



Here a is the tension of the string and £x is a coefficient 
of proportionality which is expected to be unity in (3 + 1) 
dimensions 



c 



x 



(12) 



Note that the range of qualitative applicability of the re- 
sult (fTU]) is dictated by purely geometrical reasons: the 
length R of the string should be larger than some critical 
value Rx re of order of the string width. In the case of 
unprojected action density, which is related to the heavy 
quark potential via the action sum rules [f| , one assumes 
that the coefficient a in Eq. (fTTj) equals to the physi- 
cal string tension. For the HP 1 -projected action density 
we also expect the same: although the sum rules are 
not applicable any longer, physics-wise it seems natural 
that a must be equal to the HP 1 -projected string tension. 
However, for the topological density the arguments break 
down and the scale of the width of the topological fluc- 
tuations in the vicinity of the string should be examined 
in the numerical simulations. 
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TABLE I: The parameters of our simulations. 



IV. DETAILS OF NUMERICAL SIMULATIONS 

Our measurements were performed on two sets (Ta- 
ble |T|) of statistically independent SU (2) gauge configu- 
rations generated with standard Wilson action. These 
thermalized gauge ensembles were used only to produce 
the embedded HP 1 a- model configurations, the proce- 
dure being exactly the same as in the second paper of 
Ref. [l8|. The physical scale for the obtained er-model 
data sets could be fixed by the requirement that cither 
the full string tension <7f u n or its HP 1 -projected counter- 
part cthp equals to the conventional value (440 MeV) 2 . 
Physically these two prescriptions are equivalent since 
it is known that the HP 1 -projected string tension coin- 
cides with the full one in the continuum limit. However, 
at necessarily finite (3 values a particular choice might 
be preferable for scaling checks. For us the prescription 
^/<7f u n = 440 MeV seems to be unnatural since after all 
the original gauge configurations were not used in the 
measurements. What is physically relevant in the present 
context is the HP 1 -projected string tension and we set the 
scale by fixing cthp to the above value. Note that since 
this prescription differs from the conventional one we use 
the dimensionless combinations like -R^/cthp instead of 
the usual fermi units. 

The string profiles may obviously be affected by the 
finiteness of the temporal extension T. To reduce 
these finite-time effects we used the conventional APE- 
smearing procedure [33| for spatial pieces of the trajec- 
tory C and checked independence of our results against 
variation of T. The geometrical setup described above 
makes it convenient to consider Wilson loops of only even 
temporal and spacial extents in lattice units. Further- 
more, the non-locality issue, Eq. (JSJ), forces us to disre- 
gard the smallest loops. Thus we confined ourselves to 
planar T x R Wilson contours with T, R = 4, 14 with 
understanding that T, R = 4, 6 data points must be used 
with care. 



V. TOPOLOGICAL DENSITY IN THE STRING 

The qualitative topological portrait of the confining 
string is revealed by the correlation function ([6]), shown 
in Figure [2] for our j3 = 2.60 data set. As is apparent 
from this Figure the quantity Qc is always negative in 
the vicinity of the string while approaching zero far away 
from quark-antiquark pair. Thus the confining flux tube 
suppresses the topological fluctuations. Note also that 
the suppression is maximal at the positions of quarks, 
which are sitting at the absolute minima of the topologi- 
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FIG. 2: (Color online) Correlation function (0 at quark- 
antiquark separation R^/anp = 1.47(3) measured on /3 = 2.60 
data set. Temporal extent of the Wilson loop is T^/cthp = 
1.05(2). 



FIG. 4: (Color online) The topological density correlation 
function ((6]) at the geometrical center of the string (h = r = 0) 
versus the quark separation i? for the Wilson loops of various 
temporal extension. 
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FIG. 3: (Color online) The longitudinal slices (r = 0) of the 
topological density correlation function (|6} at various quark- 
antiquark separations R at j3 = 2.60. Only the positive half 
(h ^ 0) of the axis is shown, lines are plotted to guide the 
eye. 
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FIG. 5: (Color online) Example of the longitudinal profiles at 
r = of the correlator © at string lengths R^/anp = 1.19(3) 
for (3 = 2.51 and R^/F^ = 1.26(2) for (3 = 2.60 and various 
time extensions of the Wilson loop. 



cal density profile (this observation is in agreement with 
Ref. dHHH). Moreover, the quarks look like the sources 
of Coulomb potential in the topological density, which by 
itself a rather striking result especially in view of the cor- 
responding action density profile to be discussed below. 

The typical longitudinal (r — 0) slices of the topologi- 
cal density correlation function (J6j) are plotted in Figure[3] 
for various lengths R of the confining string. The min- 
ima of the topological density are reached in the vicinities 
of the quark and anti-quark, h m i n ~ R/2. As one goes 
outward the string, h > h nimi the quantity Qc rapidly 
approaches zero indicating that topological charge den- 
sity quickly comes to its vacuum expectation value out- 
side the string. To the contrary, if one moves toward 
the string center, h — > 0, the magnitude of Qc remains 
negative slowly approaching its limiting value at h = 0. 



The degree of the suppression of the topological den- 
sity at the geometrical center of the string r = h = 
varies with the string length. This property is visualized 
for both our data sets and for various temporal exten- 
sions T of the Wilson loops on Figure [4] One can see 
that the topological density at the string center tend to 
its vacuum value as the string gets longer. This feature 
is in agreement with the picture of the fluctuating string: 
as the string gets longer its transverse fluctuations in- 
crease in amplitude and the string gets more dispersed 
in the space. As the consequence, the correlation of all 
quantities (including the topological charge density) with 
the string gets weaker at the geometrical center with in- 
creasing string length. 

Another information imprinted in Figure[3]is small sen- 
sitivity of the central value of the topological density cor- 
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FIG. 6: (Color online) The distributions of the topological 
charge fluctuations © in the transverse direction r around 
the geometrical string center h = 0. For the sake of clarity 
the data points were y-shifted, and the corresponding Qc ~ 
levels are shown by thin horizontal lines. Upper and lower 
panels refer, respectively, to the /3 = 2.51 and /3 = 2.60 data 
sets. Best fits by the Gaussian function (|1U[1 are shown by 
thick lines. The string lengths R, the temporal extensions T 
and the reduced % 2 values are indicated as well. 



relator © against the variation of the temporal extent of 
the Wilson loop T. The typical dependence of the string 
profile for h > on the time extension of the Wilson 
loop is presented on Figure [5] for the quark-antiquark 
separations Ry/aup = 1.19(3) (data set (3 = 2.51) and 
R^/W^ = 1.26(2) (/3 = 2.60). It is seen that the string 
profile is indeed robust against the time passed between 
creation and annihilation of the quark pair. Therefore 
we conclude that within our accuracy the dependence of 
the results on the time extension T is negligible. 

In order to check the validity of the Luscher for- 
mula (HHJ) for the width of the topological charge fluc- 
tuations around the string we fitted the transverse slice 
of the correlation function ([6]) at the string center h = 
by the Gaussian ansatz (jTUJ) . We find that the scaling 
law (JTUJ) well describes both our data sets as one can 
see from various examples of the fits shown in Figure [5] 



FIG. 7: (Color online) The squared string width Sq vs. the 
string length R for various time extensions of the Wilson con- 
tour. Data points had been obtained with (3 = 2.51 (top) and 
[3 = 2.60 (bottom) sets. The lines represent the best fits by 
the function (|lip (the details are given in the text). 



The only notable deviations from (JTUJ) was observed for 
small string lengths, R^/obp % 0.7, which, however, is 
to be expected since Eq. (fT0| is valid only for large R. 
Therefore in terms of the topological density the string 
formation length, at which the asymptotic behavior (jTUJ) 
settles in, could be estimated as R^/a^p m 0.7 and is 
in agreement both with the literature [l| and with the 
corresponding result coming from the action density pro- 
files (see below). In turn, the Gaussian fits allow to ob- 
tain the corresponding parameters. On the bottom panel 
of Figure [7] we plot the squared string width Sq (i?) ob- 
tained for (3 — 2.60 data set as a function of the string 
length R for various temporal extensions T of the Wil- 
son contour. One obvious thing to note here is that 
8q(R) data points are indeed compatible with logarith- 
mic dependence (fTTj) in the whole R interval for small 
T V ^5F = 0.42(1), 0.63(1), while for 7\/^hp > 0.84(2) 
they start to deviate from (fTTj) for large string lengths. 
In fact, for T^/U^ = 0.84(2), 1.05(2) the total set of R 
data points corresponds rather to the linear dependence 
of Sq on R, which is presumably a manifestation of fi- 
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nite size effects. Indeed, the considerations [l[ showing 
that finite size corrections are negligible even for largest 
loops are not applicable in our case: for HP 1 er-model 
induced gauge fields the relevant center symmetry trans- 
formations do not exist. Therefore the finite volume ef- 
fects are expected to be much more severe than they 
are in usual approaches and indeed at largest available 
T^/cthp = 1-26(2) the squared string width increases with 
R even more rapidly showing somewhat irregular behav- 
ior for large quark- antiquark separations. We attribute 
the deviation of the string width from logarithmic law 
to the boundary conditions in the temporal and spatial 
directions: the long enough string interacts with its tem- 
poral and spatial replicas. The strength of this artificial 
interaction is rising with increase of the area of the string 
worldsheet achieved at large R and T. The assumption 
of strong finite volume effects present for j3 = 2.60 con- 
figurations is justified below (see also Section |VI]) . How- 
ever, for small T the finite volume corrections are surely 
negligible, which permits us to fit Sq(R) to Eq. (fTT|) in 
the range R^/a^ > 0.7 for T^Sip = 0.42(1), 0.63(1) 
and in the interval 0.7 < -R^/ohp < 1.2 at T^/cthp = 
0.84(2), 1.05(2). The lower bound on R is due to Eq. © 
and because the Gaussian fit IjlOp is not adequate at small 
distances resulting in rather large x 2 values. The out- 
come of these fits is used below. 

In order to check finite volume effects we performed 
the same calculations on our (3 = 2.51 configurations, 
for which the physical volume is much larger. The re- 
sults are presented on the top panel of Figure from 
which it is apparent that the data points are much more 
stable in this case. Namely, the dependence Sq(R) fol- 
lows the logarithmic law (fTTj) up to the time scale of 
order T^/cthp < 1-2. However, at even larger T we 
again see notable deviations from Eq. (fTl"]) the source of 
which is not so evident for us at present. Most probably, 
the errors bars are underestimated for the corresponding 
largest Wilson loops and much more statistics in needed 
to confirm flTTJ) at T^/a^ > 1.2, R^/a^ > 1.4. 

To summarize, our numerical data on the string width 
obtained from the topological density correlation func- 
tion indeed favors the logarithmic widening (fTl]) of 
the confining string. Unfortunately j3 — 2.60 data set 
has relatively small physical volume which disturbs the 
expected dependence Sq(R). Anyway, for not so large 
time extent of the Wilson loops logarithmic string broad- 
ening is clearly seen. As far as the measurements at 
(3 = 2.51 are concerned, the corresponding physical vol- 
ume is much larger and the logarithmic widening is con- 
firmed for much larger Wilson loops. 

The fits of the Gaussian string profiles (fTTJj) allow us to 
estimate the dimensionless parameters C,q and R^ re of 
the Liischer scaling law (fTT[) . According to the theoreti- 
cal expectation the parameter must of order one pro- 
vided that the width 5q is given in terms of the physical 
string tension, while there are essentially no theoretical 
predictions for the string core width Rq tc . However, the 
relevant question for both quantities is whether they are 
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FIG. 8: (Color online) The coefficients Cq/I^chp ■ a 2 ] (top) 
and Rg Te /a (bottom) of the Liischer scaling law pip as de- 
termined from Wilson loops of various time extension T for 
both our data sets. Note that all quantities are given in the 
lattice units. 



given in physical units or are of order lattice spacing. In 
the upper and lower panels of Figure[5]we plot the param- 
eters Cq/I^hp • a 2 ] and i?g rc /a, respectively, obtained 
by the logarithmic fits (fTT[) as a functions of the temporal 
string extension T. It is apparent that T-dependence is 
indeed mild in both cases. At the same time, it is es- 
sential that the parameters are expressed in the units of 
UV cutoff, which we found the only possibility to make 
(3 = 2.51 and (3 — 2.60 data points qualitatively consis- 
tent. Hence the string width and the size of the string 
core, as they are seen by the topological fluctuations, are 
likely to be at the scale of the lattice spacing. Unfortu- 
nately, this result is mainly qualitative and is not fairly 
convincing because we consider only two lattice spacings 
while the logarithmic fits are weakly sensitive to varia- 
tions of i?g rc . For these reasons we refrain from definite 
conclusions since more data is needed to clarify the issue. 
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FIG. 9: (Color online) The action density correlation func- 
tion measured on f3 — 2.51 configurations for the string 
length Ry/aiip = 1.67(4) and for the temporal extension of 
the Wilson loop T^/^ = 1.19(3). 



FIG. 10: (Color online) The longitudinal profile (r = 0) of 
the action density correlator Sc obtained with = 2.51 data 
set. The positive half (h 0) of the axis is shown for various 
indicated string lengths R. The lines are plotted to guide eye. 



VI. ACTION DENSITY IN THE STRING 

In this Section we repeat the above investigations of 
the string profiles using the HP 1 -projected action density 



where the inessential normalization factor was omitted. 
Note that in this case the obtained statistics is much 
better allowing more precise analysis. The typical dis- 
tribution of the action density around the string is pre- 
sented on Figure [9] It is clear from the figure that the 
action density is indeed suppressed in the vicinity of the 
string in agreement with the analogous result for unpro- 
jected fields [6J] . At the same time, there are a few notice- 
able differences between the action density profile of the 
string and the corresponding topological charge correla- 
tion function, Figure O First, the positions of the test 
quark and the anti-quark are explicitly noticeable in the 
string topological profile as these positions correspond to 
the two global minima of Qq- However, the location the 
test charges in the action density correlator is not seen 
clearly. Second, the transverse distribution of the action 
is seemingly wider compared to the one of the topological 
charge. Finally, the data for the action comes with much 
less noise compared to one of the topological distribution 
at the same numerical statistics. 

In order to illustrate the above observations we plot 
in Figure [TU] the longitudinal profiles of the action den- 
sity and compare them with the same quantity for the 
topological charge plotted in Figure [3] One notices that 
the global minimum of the quantity Sc shifts with in- 
creasing distance between test color charges, being lo- 
cated at the string center point for relatively small R. 
As R enlarges the minimum shifts roughly to the string 
end-point, however, it remains much wider compared to 
that of the topological density. This property is also il- 
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FIG. 11: (Color online) The action density at the center 
of the string (h=r=Q) vs. the string length R at various time 
extensions T. 

lustrated in Figure [TT] in which the correlation function 
© at the geometrical center r = h = of the string is 
plotted as a function of the string length R. We conclude 
qualitatively that the location of the quarks becomes vis- 
ible in the longitudinal action density profile at string 
length about R^aup w 1.0. Note that the apparent dis- 
crepancy in this number between (3 = 2.51 and f3 — 2.60 
data points is presumably due to the finite volume effects 
discussed in previous section. 

As it was done in the case of the topological charge, 
we study the transverse slices of the action density corre- 
lation function ([7]) at the geometrical center h = of the 
string. Generically, one expects that the Gaussian distri- 
bution (flC)]) should adequately describe the data provided 
that the string length is large enough. It turns out that 
indeed the transverse action density distribution around 
the string obeys (JTUJ) for R^aup > 0.7 as is illustrated 
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FIG. 12: (Color online) Transverse slices of the action density 
correlation function @ at the geometrical center of the string 
h — for various temporal T and spacial R extensions of 
the Wilson contour. Upper (lower) panel refers to (3 = 2.51 
(/3 = 2.60) data set. For clarity reasons data points were y- 
shifted by —0.1, —0.2 and —0.3 (the corresponding Sc = 
levels are shown by thin horizontal lines). The fits of the 
distributions by the Gaussian function (|10|l are shown by thick 
lines. 

on Figure [T3J For smaller string lengths the data shows 
significant deviations from the Gaussian ansatz indicat- 
ing that the string formation length is not reached yet. 
Therefore, the string width could be extracted reliably 
from fits only for R-Ja^p > 0.7 and in the considera- 
tions below we restrict ourselves to these distances. 

As far as the parameters of the Gaussian distribution 
are concerned, we plot the corresponding squared width 
of the action density profile obtained at various R and 
T on Figure 1131 As is apparent from that Figure, the 
data points for 8%{R) at each j3- value confirm the loga- 
rithmic string widening for all available T in the entire 
range of quark-antiquark separations beyond the string 
formation length. Moreover, the function 5 2 S {R) appears 
to depend very mildly on T. However, basing on our ex- 
perience with the topological density distribution we ex- 
cluded from our analysis the largest T data points, which 
presumably suffer from finite volume corrections. 



FIG. 13: (Color online) The squared string width <5| vs. the 
string length R for various time extensions T of the Wilson 
contour. Data points corresponds to /3 = 2.51 (top) and /3 = 
2.60 (bottom) sets. The lines represent the best fits to {TT}. 



Next we estimate the parameters Cs and i?™ re enter- 
ing the logarithmic scaling law (fTTjl . To this end, we 
fitted the Sg(R) dependence to Eq. (fTTj) in the range 
Ry^HP > 0.7 for T values, for which the finite vol- 
ume corrections are expected to be small, namely, for 
T^/cthp < 1-0 on (3 = 2.51 configurations and for 
7\/oip = 0.42(1), 0.63(1) with (3 = 2.60 data set. For 
larger temporal extensions of the Wilson loops we re- 
stricted the fitting range from above, R^J<Jm> ^ 1.3, 
although we checked that inclusion of largest R points 
changes the results at most by ~ 10%. The outcome of 
the fits is presented on Figure [TH As might be expected, 
the amplitude of the logarithmic widening seems to be 
given indeed in physical units, as is illustrated on the top 
panel of Figure [Mj Namely, if we assume that the string 
width scale er, Eq. (fTT|) . equals to the projected string 
tension cthp then the corresponding dimensionless coeffi- 
cient £s appears to be spacing independent and is about 
unity. However, the data for string core width favors 
again the constancy of Rg ore in lattice units, the value of 
which is close to one obtained from the topological den- 
sity profiles (see the bottom panel of Figure [T4"|) . Note 
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FIG. 14: (Color online) The coefficients Cs (top) and R c s orc /a 
(bottom) of the Liischer scaling law (|ll|l as determined from 
Wilson loops of various time extensions T calculated at both 
available data sets. Note that £s is shown in physical units 
with the normalization a = cthp , while R^ la is plotted in the 
lattice units. 



that due to the generically expected small sensitivity of 
logarithmic fits to R s ore , our statistics and the spacing 
variations considered might be insufficient to determine 
accurately the width of the string core. Hence at present 
we cannot fairly insist on the relation R c g re ~ a, new 
measurements are needed to confirm it. 

To summarize, the measurements of the action density 
correlation function ([7]) confirm the logarithmic widen- 
ing of the confining string. Moreover, we were able to 
demonstrate that the string broadening is in agreement 
with general theoretical expectations and is qualitatively 
compatible with the bosonic string picture of the con- 
fining flux tube. However, the data for the string core 
width, where the string widening settles in, is less con- 
clusive and favors rather small value of Rg ore of order 



lattice spacing. In particular, even if one would insist 
that R c g re is given in physical units, then its numerical 
value is about 



r>corc 



0.1 fm. (14) 
Note that Eq. (fT4"|) is in qualitative agreement with the 



literature. Indeed, it had been shown in Refs. [34|, l35| 
that the transition from perturbative to stringy descrip- 
tion of quark-antiquark pair happens at rather small dis- 
tances of order 0.2 0.3 fm [3J]. The early start of the 
string regime was also found in Ref. [36[ in compact U(l) 
gauge theory in four space-time dimensions. It could well 
be that the early onset of the string fluctuations is a 
universal feature of the theories possessing string excita- 
tions. 



VII. DISCUSSION AND CONCLUSION 

In this paper we investigated the topological and the 
action density profiles of the confining string in pure 
SU(2) Yang-Mills theory on the lattice using the method 
of the quaternionic projective cr-model embedding. In- 
deed, the main difficulty in the conventional considera- 
tions is caused by UV dominated zero-point fluctuations, 
which give identically vanishing but numerically destruc- 
tive contribution to the correlation functions of the topo- 
logical/ action density and the confining flux tube. On the 
other hand, the HP 1 a- model embedding and correspond- 
ing gauge covariant HE^-projection procedure allows to 
filter out essentially the only leading power divergences 
from local observables in question, leaving intact the sub- 
leading and the IR dominated contributions. At very 
least this property HP 1 embedding is known to hold in 
numerical simulations and hence allows to obtain rather 
accurate results for the topological and action density 
profiles of the confining string. It is worth mentioning 
that the known disadvantages of HP 1 embedding are ex- 
pected to play no role in the problem considered. Indeed, 
the finite non-locality range of the approach should be 
inessential for large enough Wilson loop, which are rele- 
vant for the present investigation. Moreover, our results 
are in agreement with the existing literature, when the 
comparison is possible. 

Our paper presents the first quantitative investigation 
of the topological density distribution in the vicinity of 
the chromoelectric string. Although it had been known 
for some time [l3|, Q] that confining string suppresses the 
topological fluctuations, we were able to investigate the 
problem in details and to show, in particular, that the 
confining flux tube in terms of the topological density 
corresponds likely to the fluctuating string. Moreover, 
we confirmed the generically expected logarithmic string 
widening and we attempted also to estimate the charac- 
teristic string scales. Our data seems to indicate that in 
terms of the topological density the confining flux tube is 
rather thin object stretched in between quark-antiquark 
pair. Although at present we cannot fairly insist on this 
picture and more data is obviously needed to confirm the 
result, nevertheless, it is rather striking and theoretically 
challenging by itself. We hope to shed more light on the 
issue in forthcoming publications. 

As far as the action density profiles are concerned, the 
use of the projected action density allowed us to obtain 
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rather accurate results confirming the string picture. In 
this case, the string formation and its logarithmic widen- 
ing is also established. Moreover, the characteristic width 
of the flux tube is shown to be given in physical units and 
the features of its broadening seem to be compatible with 
the effective bosonic string description. 
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